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Here, we show that the conductivity of conductor-insulator composites in which electrons can 
tunnel from each conducting particle to all others may display both percolation and tunneling (i.e., 
hopping) regimes depending on few characteristics of the composite. Specifically, we find that the 
relevant parameters that give rise to one regime or the other are D/£ (where D is the size of the 
conducting particles and £ is the tunneling length) and the specific composite microstructure. For 
large values of D/£, percolation arises when the composite microstructure can be modeled as a reg- 
ular lattice that is fractionally occupied by conducting particle, while the tunneling regime is always 
obtained for equilibrium distributions of conducting particles in a continuum insulating matrix. 
As D/£ decreases the percolating behavior of the conductivity of lattice-like composites gradually 
crosses over to the tunneling-like regime characterizing particle dispersions in the continuum. For 
D/£ values lower than D/£ ~ 5 the conductivity has tunneling-like behavior independent of the 
specific microstructure of the composite. 
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I. INTRODUCTION 

The conductivity of a conductor-insulator composite 
material is characterized by a strong dependence on the 
volume fraction <fi of the conducting phase present in 
the system, and is generally understood as a percola- 
tion phenomenon arising from the electrical connectiv- 
ity between neighboring or adjacent conducting particles. 
Specifically, percolation theory considers the conducting 
particles as either electrically connected, with some fi- 
nite inter-particle conductance, or disconnected The 
introduction of this sharp cut-off implies then that be- 
low a specific fraction <fi c (the percolation threshold) the 
conductivity is zero because there is no sample-spanning 
network of connected particles, while for <p > <p c such 
network is formed and the conductivity increases as 
~ (4> — C )', where t is a critical exponent. 

It is easy to see that the notion of a sharp cut-off ap- 
plies well to composites made of large (of the order of 
one micron or more) conducting particles, for which two 
particles can be considered electrically connected only if 
they essentially touch each other, it is less clear for the 
case in which the conducting particles have sizes lim- 
ited to a few nanometers^ In that situation even if the 
particles do not physically touch each other, their mean 
separation in the composite is such that electrons can 
still flow via tunneling processes from one particle to the 
other. The resulting tunneling conductance decays ex- 
ponentially with the inter-particle distance over a char- 
acteristic tunneling length £ which is of the order of a 
fraction to a few nanometers depending on the material 
properties. Since the tunneling decay does not imply any 
sharp cut-off, the basic assumption of percolation theory 
in describing composites made of nanometric conducting 
fillers is not justified a priori. In fact, it turns out that 
nanocomposites in the dielectric regime, for which the 
conducting particles do not touch each other, arc better 
explained by inter-particle tunneling, with no imposed 



sharp connectivity criterion, than by the classical per- 
colation theoryi^i The resulting conductivities of these 
systems follow therefore hopping-like (or tunneling-like) 
behaviors,— ~— with no critical percolation thresholds as a 
function of 4>. 

The possibility of having pcrcolation-likc or tunneling- 
like regimes depending on the size D of the conducting 
particles compared to the tunneling length £, as sug- 
gested above and in Rcf. does not appear to have 
been further elaborated in the literature, despite of its 
relevance to the understanding of transport properties in 
conductor-insulator composites. In the present paper we 
address this issue by considering a global tunneling net- 
work (GTN) model of conductor-insulator composites, - 
where each conducting particle is connected to all others 
via tunneling processes. We show that this model permits 
to treat the percolation and the tunneling regimes on 
equal footing, and that one can switch from one regime to 
the other depending not only on £)/£ but, most notably, 
also on the specific distribution of the conducting phase 
in the composite. We illustrate this behavior by con- 
sidering two idealized realizations of conductor-insulator 
composites: a lattice model, where conducting spherical 
particles of diameter D occupy randomly a fraction of the 
sites of a regular lattice, and a continuum model, where 
the conducting particles are dispersed with an equilib- 
rium distribution in a continuous insulating medium. By 
using both the effective medium approximation (EMA) 
and Monte Carlo (MC) calculations, we find that the con- 
ductivity of the continuum model has a tunneling-like be- 
havior independent of the value of D/£, while the lattice 
model displays a percolation behavior of the conductiv- 
ity only for sufficiently large D/t; values, which gradually 
crosses over to a tunneling-like regime as D/£ decreases. 
For £)/£ < 5 we show that the conductivity of the lat- 
tice model is basically indistinguishable from that of the 
continuum case, establishing therefore a crossover point 
from percolation to tunneling behaviors. 
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This paper is organized as follows. In Sec.|TI]we intro- 
duce the GTN model, we formulate the EMA for both 
lattice and continuum models and we calculate the result- 
ing EMA conductances. In Sec. IIIII we present our MC 
calculations for both the lattice and continuum models, 
and in Sec. IIVI we discuss the crossover between percola- 
tion and tunneling regimes. Discussions and conclusions 
are given in Sec. [V] 



II. MODEL AND EFFECTIVE MEDIUM 
APPROXIMATION 

The GTN model for conductor-insulator composites is 
defined by considering n identical conducting particles 
contained within a volume V . This defines the volume 
fraction <p = pv of the conducting phase, where p = n/V 
is the particle density and v is the volume of a single 
particle. For the purpose of the present work we limit the 
analysis to the relatively simple case of spherical particles 
of identical diameter D, so that v = irD 3 /6. Next, we 
assume that the tunneling conductance between any two 
particles centered at and Vj is given by: 



g(n 



go cxp 



(i) 



where go is a constant "contact" conductance which in 
the following we shall set equal to unity, £ is the char- 
acteristic tunneling length, and =| — Vj | is the 
distance between two sphere centers. For simplicity we 
further assume that selective tunneling mechanisms aris- 
ing from excitation energies can be safely neglected (as 
in the case for nearest neighbor hopping^). 

The set of all tunneling conductances of Eq. (JlJ defines 
a resistor network whose conductivity depends on £ and 
-D, as well as on the volume fraction <j) and on the specific 
distribution of the particle centers. As shown below, it 
turns out that all these dependencies are well captured by 
the (single bond) EMA applied to the tunneling resistor 
network^ This is given by the solution of the following 
equation: 



E 
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g(nj) + [z l z j /(z z + Zj) 



= 0, 



(2) 



where g is the effective bond conductance and i and j run 
over the positions of the n particles. In the above expres- 
sion, the symbol (. . .) denotes the statistical average over 
all realizations of the n-particle system. Furthermore, Zi 
and Zj are the coordination numbers of the two ends of 
a resistor g(rij), i.e., given a particle at i (j), Zi (Zj) 
is the number of particles that are electrically connected 
to i (j). Since it is assumed that all particles are con- 
nected to each other regardless of their relative distances, 
Zi = Zi = 



1. In this way, Eq. ([2]) reduces to 



E 



g(i 



S (ry) + [(n-l)/2-l]g 



= 0. 



(3) 



By noticing that J2i^j = n ( n ~ 1) an d by using Eq. ([l}, 
we can write Eq. @ more conveniently as 



n^.g* exp[2(r y -£>)/£] + !_ 



= 2, 



(4) 



where for large n the dimensionless conductance g* = 
[(n — l)/2 — l]g/go — ng/2go coincides with the conduc- 
tance between two random nodes of the network with all 
bond conductances equal to g/go^ 

By multiplying each term of the summation over i,j 
by fdr5(r — r^ ) = 1, where the integration is extended 
over the entire volume V, the left-hand side of Eq. d3| 
becomes 



1 



ff*exp[2(r y -D)/£] + l 



= /^(^E^-^)) ff . exp[2(r l 



D)/t] + r 



(5) 



The above expression can be rewritten in terms of the 
pair distribution function gi (r) by noticing that, by 
definition^ 



P52 



so that Eq. ((4]) reduces to 

4irr 2 pg 2 {r) 



dr 



g* exp[2(r -£>)/£] + ! 



= 2, 



(6) 



(7) 



where g 2 (r) = Jdfl g 2 (r) / 4ir is the radial distribution 
function (rdf) and the upper limit of the integration has 
been set to infinity because of the exponential decay of 
the integrand. Note that Eq. [7] is similar to the result 
obtained in Ref. [TH. 

An alternative and useful version of Eq. ((7]) can be 
obtained by introducing the function 



W(r) = 



with 



f 



g* exp[2(r - £>)/£] + 1 exp[2(r - r*)/£] + 1 ' 

(8) 



r* =D + -]n( — 
2 W 



and integrating Eq. ([7]) by parts: 

dW{r) 



drZ(r) 



dr 



where 



Z[r) = / dr'4irr' 2 pg 2 (r'), 



(9) 



(10) 



(11) 



is the cumulative coordination number function (i.e., the 
number of spheres whose centers arc within a distance r 
from the center of a given sphere). 
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FIG. 1: (Color online) (a) Radial distribution function gi{v) 
for the lattice (dotted lines) and the continuum (solid lines) 
models. For the lattice model the vertical dotted lines in- 
dicate the positions of the delta-peaks centered at Rk for 
k = 1, 2, 3, . . .. The g2{r) curves for different <j) values of the 
continuum model have been obtained by using the results of 
Ref . [Til, (b) Cumulative coordination number Z(r), Eq. {TTJ, 
for the lattice (dotted lines) and the continuum (solid lines) 
models for different d> values. 



in Eq. © reduces to 

R^O 

which leads to an rdf of the form: 



W2(r) = |- j£ 6 ( r - R k), 

fe=l,2,... k 



(12) 



(13) 



where N k is the number of the fc-th nearest neighbors 
being at distance Rk from a reference particle set at the 
origin. The characteristic feature of this fractionally oc- 
cupied lattice model is therefore that its rdf, shown in 
Fig. QJa) by dotted lines, is given by delta-peaks whose 
positions do not change with p. As we shall see in the 
following, this feature is directly related to the appear- 
ance of a percolation behavior of the EMA conductance 
g* for large D/£ values. 

The second model of particle distribution is given by 
an equilibrated dispersion of impenetrable spheres in the 
continuum. Typical examples of the resulting 172 (?") for 
this case, obtained by using the formula provided in 
Ref. [3 are shown in Fig. [2a) by solid lines and for dif- 
ferent values of the volume fraction 0. As opposed to the 
lattice case, the rdf of an equilibrium dispersion of hard 
spheres is continuous and nonzero in the whole r > D 
range, which we anticipate here to be the pre-requisite 
for an hopping behavior of the composite conductivity. 

In Fig. HJb) we compare the coordination number func- 
tion Z(r), Eq. (fTT|) , for the lattice (dotted lines) and con- 
tinuum (solid lines) models. Contrary to the continuum 
model in which Z(r) increases smoothly with r, the lat- 
tice Z(r) has a step-like increase. At large r however the 
main contribution of the lattice Z(r) goes as (4/3)irpr 3 , 
meaning that the distribution of spheres in a sparsely 
occupied lattice is basically that of point particles in the 
continuum. 



A. Lattice and continuum models 

As it is clear from the EMA equations ([7]) and (fTCTj) , all 
the informations on the spatial distribution of the con- 
ducting particles is contained in the rdf <?2 (f), whose de- 
pendencies on r and p govern the behavior of the EMA 
conductance g*. Here, we consider two possible real- 
izations of a conductor-insulator composite which cor- 
responds to two rather extreme forms of gi (r). 

In the first case, we consider a simple cubic lattice with 
a lattice constant equal to the sphere diameter D, where 
only a fraction p of the lattice sites is occupied, randomly, 
by the spherical conducting particles, while the remaining 
fraction 1 — p is occupied by insulating spheres of equal 
diameter (as in the Scher and Zallen model of Ref. Il3f ) . 
The i-th conducting particle occupies therefore the po- 
sition r.; = R with probability p = n/N, where R is 
the direct lattice vector running over all N sites of the 
cube. In this way, the pair distribution function defined 



B. EMA conductance 

Let us start by considering the solution of the EMA 
equation for the fractionally occupied cubic lattice case. 
By using Eq. ([L3"]) we find that Eq. ([7]) reduces to: 

„ V — = 2 (14) 

fe =^...3* ex P[ 2 (^-^] + l ' 1 ' 

In the limiting case of very large particle sizes such that 
D /£ — > 00, only the first (fc = 1) nearest neighbors with 
Ri = D and N\ = 6 contribute to the summation in 
Eq. (fTl")) . which reduces to 6p/(g* + 1) = 2. Hence, in 
this limit, the resulting conductivity g* = 3(p — 1/3) 
has a percolation behavior and vanishes at the critical 
fractional = 1/3. For finite D/£ values, the next nearest 
neighbors also contribute to the network conductivity. 
For example, by assuming that D/£ is large enough to 
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retain only the terms up to the second nearest neighbors 
(N 2 = 12 and R 2 = V^D) we get that Eq. flU) becomes 



6p 



12p 



g* + l g*exp[2D(V2- !)/£] + ! 



= 2, 



(15) 



whose solution still behaves as g* ~ p — p\ for p > p\ 
but remains finite, albeit exponentially small, at lower p 
values. This is because the tunneling to the 2-nd near- 
est neighbors vanishes only at p = pa = 1/9. Con- 
sidering the whole set of neighbors when p decreases, 
one finds that the solution of Eq. (|T4|) becomes char- 
acterized by a monotonous decrease of g* punctuated 
by sharp (for small D/£ values) drops at p = pi~ with 

Pk = 2/ {J2k'=i Nk')- Furthermore, in the vicinity of each 
Pk , the conductance follows the EMA power law behavior 
— (P-Pfe) for p > p k . 

This feature is illustrated in Fig. (5Ja) where we show 
numerical solutions of Eq. flT4"|) for different values of 
D/£. For D/£ = 100, the first percolation transition 
at pi = 1/3 ((f>i = pi7r/6 ~ 0.175) causes a drop of 
the conductivity of about 36 orders of magnitude (not 
shown in the figure) compared to the conductivity at 
the close packing fraction tt/6 that corresponds to the 
lattice sites that are all occupied by the conducting 
spheres. Such drop of the conductivity is well beyond 
the maximum range of conductivity values measured in 
real composites^ and so, for all practical purposes, the 
D/£ = 100 case behaves as (0— 0i)* with t = 1 being the 
EMA critical exponent. For £)/£ = 25 the first transition 
at <f>x leads to a drop of g* of only 8 orders of magnitude, 
and the transitions at lower values of <f> are clearly vis- 
ible in the figure. For lower values of D/£, the drops 
of g* are further reduced and the transitions are much 
smoother, due of course to the fact that for these D/£ 
values the probability of tunneling to neighbors that are 
farther apart is enhanced. For D/£ = 5 the variation of 



g* is so smooth in the whole range of <f> that the under- 
lying lattice structure can be considered as completely 
irrelevant. 

In contrast with the lattice case, the EMA conductance 
g* resulting from an equilibrium distribution of conduct- 
ing spheres in the continuum does not display any perco- 
lation behavior even for large D/£ values. This is shown 
in Fig. [5Jb) where solutions of Eq. J7]), with g<i (r) as 
given in Ref . [HI, are plotted for the same D/£ values of 
Fig. Ha). The lack of percolation behavior in this case 
is due to the fact that the corresponding rdf [solid lines 
in Fig. [TJa)] is always nonzero for r > D and it does not 
vary much even for <j) values close to the packing frac- 
tion of the simple cubic lattice (cj) = tt/6 ~ 0.524). To 
see how the corresponding tunneling behavior arises, let 
us consider the EMA in the form of Eq. (|10p. Since, 
as shown in Fig.QJb), the coordination number function 
Z(r) is a smooth increasing function of r/D and, given 
that —dW/dr is peaked at r = r* with spread ~ we 
can approximate Eq. (|10[) by 



Z(r*) = 2, 



(16) 



so that, by applying Eq. (j9)), the EMA conductance be- 
comes 



g = exp 



2(r* - D) 
I 



(17) 



where r* is such that Eq. (|16|) is satisfied. In passing, 
we note that Eqs. ([To) and (fT?) represent the EMA 
equivalent of the critical path approximation (CPA) of 
Refs. ||||9pi For sufficiently large r* (i.e., small <?*) the 
coordination number goes as Z(r*) ~ (4/3)irp(r*) 3 and 
from Eq. (fl7|) we obtain that 



exp 



2D / 0.63 



1 



(18) 
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which coincides, if the coefficient 0.63 in the exponent 
is replaced by 0.7, with the low-density hopping behav- 
ior as obtained from the CPA.— In the same low-density 
limit, Eq. (fTS"]) (with a slightly different coefficient in the 
exponent) has been derived also in Ref. [lU by using an 
EMA-based procedure similar to the one presented here. 

The hopping-like (or tunneling-like) dependence of g* 
for the continuum model is further illustrated in Fig.[2Jc) 
where the results of Fig. [2jb) have been re-plotted as a 
function of (dotted lines). In Fig. [JJc) we show 

for comparison also the g* versus curves (solid 

lines) of the lattice model results of Fig. [^a), which 
tend, with the decrease of D/£, to the tunneling behav- 
ior of the continuum model. An asymptotic equivalence 
between the two models is indeed expected because, as 
noted above, the underlying microstructures become ir- 
relevant for small £>/£■ Mathematically, this can be seen 
from Eq. ([TU| since —dW/dr averages Z(r) within a dis- 
tance ~ £ around r* , and thus the differences between 
the lattice and the continuum models are blurred for suf- 
ficiently large £. Note also that for the lattice model 
Z(r*) ~ (4/3)-7rp(r*) 3 for large r* , implying that g* tends 
asymptotically to the results given by Eq. (|T8j) , as shown 
by the solid lines in Fig. [2jc) . 



III. MONTE CARLO RESULTS 

Let us turn now to a comparison of our EMA results 
of the previous section with full MC calculations of the 
conductivity for both the lattice and continuum mod- 
els. For the continuum we have generated equilibrium 
distributions of impenetrable spheres inside a cubic cell 
with periodic boundary conditions followingthe proce- 
dure that wc outlined previously in Rcfs. |j|5|. For the 
lattice model, we simply populated a given fraction of 
the cubic lattice sites with conducting spheres with di- 
ameter equal to the cubic lattice constant. 

To calculate the conductivity resulting from the lat- 
tice and continuum models, we ascribed to each pair of 
particles the tunneling conductance of Eq. ([1} , and per- 
formed a numerical decimation of the resulting resistor 
network4i5iiS To reduce computational times of the dec- 
imation procedure we introduced an artificial maximum 
distance between the particles beyond which the result- 
ing bond conductance can be safely neglected4^ 

The results of the calculated conductivity are plotted 
in Fig. [3] for the same parameter values that were used 
for Fig. [2j Each symbol is the outcome of Nr = 50 
realizations of systems of Np ~ 1000 spheres. The loga- 
rithm average of the results was considered since, due to 
the exponential dependence of Eq. (TTJ, the distribution 
of the computed conductivities was approximately of the 
log-normal form. 

It is interesting to notice that there is an overall quasi- 
quantitative agreement between the EMA and the MC 
results, for both lattice and continuum models, mean- 
ing that the EMA formulation of Sec. [TT| captures well 



the physics of the problem. For example, the first per- 
colation threshold at p\ = 1/3 (0i = npi/6 — 0.175) 
obtained from the EMA on the cubic lattice is very close 
to the critical value p\ ~ 0.3116 (</>i ~ 0.163) for the site 
percolation problem on the simple cubic lattice*^ 2 -*^ An 
important expected difference is however found in the re- 
gion cj> > 4>i where the MC conductivity for the lattice 
case should follow a (4> — </>i)' dependence with a critical 
exponent t ~ 2 instead of the EMA exponent of t = l. 2 
This is indeed verified in Fig. 2] where the MC conduc- 
tivities for D/£ = 25 (a) and D/£ = 15 (b) are both 
fitted with t ~ 1.56 ± 0.04, which is slightly lower than 
t ~ 2 due to finite size effects (for the same reason we 
find (f>i ~ 0.18 instead of 4>i ~ 0.163). Similar critical 
behaviors are expected in the vicinity of all consecutive 
percolation thresholds 4>k with k > 1, but due to the lim- 
ited number of particle densities considered in our MC 
calculations and the close proximity of successive pk val- 
ues we have been able to fit only the D/£ = 25 case in 
the vicinity of the second percolation threshold. Wc have 
found t = 1.76 ± 0.15, which agrees within errors with 
the exponent at </>i, and <j>2 — 0.074 (p 2 — 0.141) which 
is close to the expected value </> 2 ~ 0.072 (p 2 — 0.137)^S 
These same values of t and <p2 appear to reproduce quite 
well the percolation behavior around the second percola- 
tion threshold also for the D/£ = 15 case [Fig. HJb)]. 

Another important agreement between the EMA and 
the MC results is the tendency of the lattice model 
conductivity to become tunneling-like as D/£ decreases. 
This is illustrated in Fig. [3jc) where the MC conductiv- 
ities of Fig. (3fa) and Fig. 0(b) are plotted as functions 
of (/) -1 / 3 . The asymptotic regime for small a of the con- 
tinuum model (open symbols) follow Eq. (fT5|) with the 
coefficient 0.63 replaced by 0.7, as obtained from the crit- 
ical path approximation applied to a dilute system of im- 
penetrable spheres^ In complete analogy with the EMA 
results, the conductivities of the lattice model [filled sym- 
bols in Fig. 113(c)] steadily tend to this tunneling regime 
as D/£ decreases, and for D/£ = 5 the MC results of the 
two models are practically indistinguishable for all values 
of 0. 



IV. PERCOLATION-TO-HOPPING 
CROSSOVER 

Both EMA and MC results point toward a substantial 
equivalence between the cubic lattice and the continuum 
models for sufficiently small values of D/£. Given that 
the two models considered represent two extreme limits 
of how the conducting particles may be ideally arranged 
in an insulator-conductor composite, such equivalence is 
important for the understanding of the transport prop- 
erties in real composites, whose microstructure is neither 
completely lattice-like nor exactly an equilibrium distri- 
bution in the continuum. We find it therefore useful to 
define a measure for the deviation between the lattice 
and continuum models in order to follow quantitatively 
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FIG. 3: (Color online) Monte Carlo conductance for (a) the cubic lattice model and (b) the continuum model for different 
values of D/£. (c) The same results of (a) and (b) plotted as a function of 0~ 1//3 . 
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FIG. 4: (Color online) The lattice Monte Carlo conductivity 
as a function of for (a) D/£ = 25 and (b) = 15. The 
solid lines are the best fits with the percolation formula a — 
aoU - 0c)*- 



how these two extremes approach each other as D/£ de- 
creases. 

To this end, we introduce the following quantity: 



lO! 



/ CTiatt. 



\ f^cont. 



(19) 



where <J\ a tt. and cr con t. are the conductivities for the lat- 
tice and the continuum models, respectively, and the 
log is the logarithm with base 10. Hence, according to 
Eq. (fTIj|) , A = M if for a certain value <7i a tt. differs 
from cr con t. by M orders of magnitude. In Fig. [5] we plot 
the maximum value A max and the mean value A mcan 
of A as calculated over the entire range of <f> considered 
and for several values of D/£. The open symbols are 



calculations while the filled symbols refer to the MC re- 
sults. As clearly seen in the figure, when £>/£ decreases 
both EMA and MC data display a steady decrease of 
A max and A mcan . For D/£ < 10, A max becomes less 
than unity, which means that eri a tt. and cr con t. differ at 
most by less than one order of magnitude in the whole 
range of </> values. We also note that for £)/£ < 10 the 
MC results display a somewhat stronger decrease than 
the EMA results. This feature does not seems to be due 
to finite size effects in the MC calculations, and we at- 
tribute it to a real deviation from the EMA results. We 
can then identify the value of D/£ ~ 5 as a crossover be- 
tween percolation and hopping regimes below which the 
conductivities of the lattice and continuum models differ 
significantly less than one order of magnitude. 

As repeatedly stressed above, the lattice and contin- 
uum models are rather idealized representations of the 
true microstructurc of real conductor-insulator compos- 
ites. However they also define two extreme boundaries 
which delimit somehow all the possible configurations 
that can be found in isotropic and homogeneous com- 
posites. For this reason, the results of Fig. [5] can be 
considered as upper boundaries of more realistic com- 
posites, and our D/£ ~ 5 result probably underestimates 
the crossover point between percolation and hopping for 
composites whose microstructurc deviates from an ideal 
lattice model. 



V. DISCUSSION AND CONCLUSIONS 

The results presented in the previous sections consis- 
tently show that the GTN model of conductor-insulator 
composites, where each particle is connected via tunnel- 
ing to all others, is capable of explaining the appear- 
ance of both percolation and hopping regimes depend- 
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FIG. 5: (Color online) Maximum value A max and mean value 
Amcan extracted from the whole <f> dependence of Eq. (|19|l and 
for several values of D/£. The open symbols are the EMA 
results, while the filled symbols refer to the MC calculations. 



ing on few characteristics of the composite. In particu- 
lar, we have identified the ratio D/t; and the composite 
microstructure as the relevant variables that control the 
switch from one regime to the other. For large values of 
percolation arises in composites whose microstruc- 
ture can approach a regular lattice that is fractionally oc- 
cupied by conducting particle, while a hopping regime of 
the conductivity is always obtained for equilibrium distri- 
butions of conducting particles in the continuum. As D/t; 
decreases, the conductivity of the lattice-like composites 
gradually loses its percolating character and approaches 
the hopping regime which characterizes the continuum 
limit. For D/£ < 5 the composite conductivity displays 
hopping behavior independent of the specific microstruc- 
ture and is practically indistinguishable from that arising 
from an equilibrium distribution of particles in the con- 
tinuum. 

Given that the values of the tunneling characteristic 
length £ do not exceed a few nanometers, our analysis 
predicts that composites whose conducting fillers have 
nanometric sizes should always display a hopping-like (i. 
e., percolation- less) behavior independent of the partic- 
ular distribution of the conducting particles in the com- 
posite. Instead, composites with conducting filler sizes 
larger than a fraction of microns should display percola- 
tion or hopping behaviors depending on the whether the 
microstructure is more lattice-like or more continuum- 
like, respectively. In view of the above, it appears then 
reasonable that dilute filler polymer-based composites 
will display hopping-like behavior also for filler diame- 
ters substantially larger than £ (i. e., in the order of some 
hundreds of nanometers) as we have verified recently in 
Ref. i. At the same time, our theory also explains the 
conductivity behavior of composites made of mixtures of 



hard conducting and insulating particles. In particular, 
we note that in dense ensembles of conducting spheres 
the arrangement in continuum systems is similar to that 
of lattices i 13 i 21 ' 22 In this respect, the results of Ref. 
on the conductivity of co-sputtered Ni-Si02 cermets rep- 
resent a nice example of tunneling-driven percolating be- 
havior on a lattice-like microstructure. Indeed, the data 
of Ref. [23] display multiple percolation thresholds as the 
concentration of conducting Ni grains is reduced, very 
much like the behavior shown in Fig. 2] for the MC con- 
ductivity of a fractionally occupied lattice model. The 
difference of about 3 orders of magnitude between the 
conductivity at the largest Ni concentration and that at 
the first percolation threshold can be reproduced by our 
model by setting D/£ ~ 10 which, by using the mea- 
sured mean Ni grain size D sa 10 nmjS leads to £ 1 
nm, which is of the expected order of magnitude (see e.g. 
Refs.ii) and compares well with the £ values extracted 
from other composites with spherical fillersi^£ It should 
be noted however that the lattice model applies only par- 
tially to the Ni-Si02 data because the nonuniversal value 
of the conductivity exponent t observed in the vicinity 
of the second percolation threshold in Ref. HH cannot be 
reproduced by our lattice MC results (see Fig. [4j. As 
to be discussed elsewhere, a nonuniversal value of t (in 
the sense that we have specified in Refs. [l8ll24l ) could be 
obtained within the GTN model by allowing a finite dis- 
persion in the distances of the second nearest neighbors. 



Before concluding it is worth to point out some lim- 
itations of the theory presented above. First, for sim- 
plicity, we have considered composites whose conducting 
fillers are given by monodispersed spheres. Even main- 
taining that for some classes of composites the shape of 
the fillers can be approximated by a sphere, the collec- 
tion of such spheres in a real composite has however some 
degree of polydispcrsivity, which can be large enough to 
make the fractionally occupied periodic lattice model an 
inappropriate description of the possible configuration of 
the microstructure. However, as stressed above, the lat- 
tice model should be regarded as an extreme deviation 
from a continuum dispersion of particles (monodispersed 
or not), so that the crossover studied in Sec. |IV] repre- 
sents an upper limit also for the case of polydispersed 
spherical particles. The second limitation is associated 
with the fact the we have limited ourselves to the case 
of isotropic particle fillers, and we have not attempted to 
analyze how percolation and hopping behaviors arise for 
anisotropic filler shapes such as rodlike or platelike parti- 
cles. An analysis of these cases would be interesting also 
in relation to the appearance of nematic phases for large 
volume fractions of anisotropic fillers and to lattice-like 
arrangements driven by attractive forces in dispersions of 
rodlike particles. 
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